Many-body spin interactions and the ground state of a dense Rydberg lattice gas 
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We study a one-dimensional atomic lattice gas in which Rydberg atoms are excited by a laser and 
whose external dynamics is frozen. We identify a parameter regime in which the Hamiltonian is 
well-approximated by a spin Hamiltonian with quasi-local many-body interactions which possesses 
an exact analytic ground state solution. This state is a superposition of all states of the system 
that are compatible with an interaction induced constraint weighted by a fugacity. We perform a 
detailed analysis of this state which exhibits a cross-over between a paramagnetic phase with short- 
ranged correlations and a crystal. This study also leads us to a class of spin models with many-body 
interactions that permit an analytic ground state solution. 
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Strong interactions competing with non-commuting 
single particle terms in a many-body quantum Hamil- 
tonian often lead to non-classical ground states. Only 
in exceptional cases analytic or approximately analytic 
results can be found. Paradigm examples are the one- 
dimensional xy-mode\ in a transverse field [T] or the cele- 
brated Affleck-Kennedy-Lieb-Tasaki spin model [2] , both 
of which have proven indispensable for the understanding 
of many-body phenomena in magnetic compounds and 
valence bond solids. Finding models of experimentally 
realizable many-body Hamiltonians with exact solutions 
is hence of fundamental interest. 

Models of many-body quantum systems with origin in 
condensed matter physics are currently implemented and 
studied in ultra cold atomic systems with great success 
[3] . Most experiments so far are carried out with ground 
state atoms, but very recent efforts exploit the unique 
properties of atoms excited to Rydberg states. These 
states offer strong and long-ranged interatomic interac- 
tions in conjunction with an extraordinarily long lifetime 
[4j [5] . This enabled remarkable experiments which stud- 
ied the coherence properties in strongly interacting Ryd- 
berg gases [51 [7] and eventually demonstrated the feasibil- 
ity to process quantum information with Rydberg atoms 
EH EH • Rydberg gases are moreover an almost ideal ex- 
perimental implementation of interacting spin systems 
such as the aforementioned xy- model |10j . This stimu- 
lated a plethora of theoretical studies investigating the 
real time evolution [TTJ Q2] as well as ground states of 
these spin models [13 15 . The latter, predominantly nu- 
merical work, revealed a variety of interesting quantum 
phases and studied the creation |T3l [14] as well as the 
melting dynamics [TS] of dynamically created crystals. 

In this work we provide an analytic study of the non- 
trivial entangled many-body ground state of a strongly 
interacting one-dimensional Rydberg gas. The strong 
interaction among excited atoms gives rise to an effec- 
tive Hamiltonian with a quasi-local three-body interac- 
tion that effectuates a set of non-commuting local con- 
straints. For certain values of the experimental param- 



eters this Hamiltonian is accurately approximated by a 
spin Hamiltonian which has an exact analytical ground 
state solution. We show that this is a consequence of the 
fact that the Hamiltonian possesses a manifold of approx- 
imate Rokhsar-Kivelson points jTHj where is assumes a 
so-called Stochastic Matrix Form [XT]. The ground state 
is a coherent superposition of all states compatible with 
all the local constraints weighted by an effective fugacity. 
We analytically explore the properties of this state which 
shows a cross-over between a paramagnetic phase with 
short-ranged correlations and a crystalline ordered state. 
Our study highlights a new perspective for creating and 
studying non-classical and entangled states with ultra- 
cold Rydberg gases. It also leads to a class of spin mod- 
els with many-body interaction whose non-trivial ground 
state solution can be found analytically. 

Our system consists of a deep one-dimensional lattice 
with L sites evenly spaced at a distance a. For con- 
venience we consider periodic boundary conditions, but 
this is no necessary requirement. Each site is occupied 
by a single atom which we treat within the two-level ap- 
proximation where every atom forms a (pseudo)spin 1/2 
particle. The atomic ground state \g) = \\) is coupled to 
a Rydberg state |r) = |f) by a laser with Rabi frequency 
n and detuning A. Atoms in Rydberg states interact via 
a power law interaction with (inverse) exponent 7 > 0. 
The Hamiltonian of this system is (within the rotating 
wave approximation for the atom-laser coupling) given 
by 

tfo = ^ + A£n fe + y£^. (1) 

k k rn>k 1 1 

Here V is the interaction strength, cr£ is a Pauli matrix 
and 7ifc = crjL<7_ is the Rydberg number operator with u\ 
being the raising and lowering operators of the fc-th spin. 
This Hamiltonian is an Ising-model with long-range in- 
teractions in a transverse and a longitudinal field. Most 
experiments up to date use atomic states that interact 
via the van-der-Waals interaction, i.e. 7 = 6. We will fo- 
cus here on this case but our approach also works for the 
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dipole-dipole interacting case with 7 = 3. The Hamilto- 
nian ([I]) was employed successfully to describe the exci- 
tation dynamics of driven Rydberg gases and has been 
proven to reflect very accurately the actual experimental 
situation [51 ITS]. 

The interaction between excited atoms strongly affects 
the excitation dynamics of the system through a mech- 
anism which is called the Rydberg blockade [19 : The 
large interaction induced energy shift makes it virtually 
impossible to excite two nearby atoms simultaneously to 
the Rydberg state, i.e. the presence of one excited atom 
blocks the excitation of atoms in its vicinity (we assume 
V 3> |A| and il V throughout). In what follows we 
will make this blockade effect manifest, which will create 
an effective three-body spin interaction in the Hamilto- 
nian. Owed to the power law decay the strongest in- 
teraction takes place between nearest neighbors and we 
assume that a strict blockade is only present between 
them. We transform the Hamiltonian ([I]) into an inter- 
action picture with respect to the nearest neighbor in- 
teraction by applying the unitary transformation U — 
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( 1 ) is the 



The first term of eq. 
only one that does not commute with U and one~bbtains 
W**U = [P fe -i + n k _ ie uv ] a\ [P k+1 + n k+1 e uv ] +h.c. 
where P k = l — n k . Note that both n k and P k are projec- 
tion operators. Since V ^> f2 one can neglect the terms 
with rapidly oscillating phases which is essentially a ro- 
tating wave approximation. We then arrive at our work- 
ing Hamiltonian 
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where the first term is the blockade-induced three-body 
interaction: The excitation of an atom to a Rydberg 
state on site k can only take place provided that both 
projectors P k -\ and P k +\ yield a non-zero value. This 
imposes a constraint such that the Hilbert space splits 
into uncoupled blocks each of which is characterized by 
the number of pairs of neighboring excited atoms. We 
will be concerned with the subspace in which there is no 
simultaneous excitation of nearest neighbors (referred to 
as physical subspace). 

In the following we will show that for a certain set 
of parameters (fl, A, V) Hamiltonian ^ possesses an 
approximate analytical solution. The decisive idea is to 
add the term H$ = J2 k Pk-iPk+A^nk + £(1 - n k )} to 
Hamiltonian ^ where £ is a real and positive parame- 
ter and subsequently subtract it. Obviously, adding and 
subtracting does not change H, but regrouping all 
terms conveniently allows us to rewrite the Hamiltonian 
as H = Eq + i/3body + H' where now each term depends 
on £. Here E — — QLt; will turn out to be the approxi- 
mate ground state energy, i?3body is a spin Hamiltonian 
with three-body interactions that has an analytic ground 



state solution and H' is a perturbation. -ff3body is given 
by 

L L 

H 3hody = jve- 1 + eE hk = n E h l h >°> ( 3 ) 



with 



Pk-iPk+x [4 + r V + e(i - n k )] (4) 



The term of h k which contains a k is proportional to the 
three-body interaction term in Hamiltonian ([2]). The re- 
maining ones are chosen such that h k becomes a self- 
adjoint operator with positive-semidefinite spectrum and 
only one non-zero eigenvalue: \J £ _1 + £. These are the 
terms that were introduced by H^. The part that is can- 
celing them is contained in the perturbation H' which 
we discuss later. At first, we construct the ground state 
|£) of the Hamiltonian i?3body This state is annihilated 
by all h k and has hence energy zero. Note, that oper- 
ators acting on neighboring sites do not commute, i.e. 
[h k , h k ±i] 7^ 0. It is thus not possible to use a local zero 
eigenstate of each h k and then form the total ground 
state by a product of these states. Instead one finds that 
the state of the physical subspace that is annihilated by 
all h k is given by 



fl (i - zp k - 1( T k + p k+1 ) \u ... i) ■ (5) 
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where is a normalization constant. This state is a 
coherent superposition of all states that have no near- 
est neighbor excitations. The probability of each state is 
weighted by the factor (£ 2 )™, where n is the total num- 
ber of Rydberg excitations in this state. This state is 
highly non-classical as it is a coherent superposition of 
all states from the physical subspace and cannot (except 
for £ = 0) be written as a product state. The existence of 
this ground state is due to the special projector property 
of each term of Hamiltonian ^ which is also known as 
Stochastic Matrix Form [17] . 

In order to calculate the normalization constant 
one has to count the number of all allowed arrangement 
of excited atoms on the lattice and sum them using the 
weights (£ 2 )™. Since there is strict nearest neighbor ex- 
clusion this sum is equivalent to the partition function 
of a lattice gas of hard-core dimers, i.e. hard objects 
that occupy two neighboring lattice sites. In the limit 
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L » 1 we obtain Z i = (1/2)(1 + a/1 + 4£ 2 ) 

that we can identify £ 2 as a fugacity. The fugacity sup- 
presses/enhances the weight of a state with n excited 
atoms or dimers by (£ 2 ) n [2Q]- We emphasize that the 
correspondence between the quantum problem and the 
dimer gas is solely formal. One striking difference is the 
range of the interaction: for the classical system only 
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nearest neighbors interact while in the quantum system 
also interaction among next-nearest neighbors occur. 

The aim is now to find a set of parameters (O, A, V) 
or a whole manifold of them such that H' is negligible 
compared to H^ ot \ y . In this case the Hamiltonian of the 
Rydberg gas is be very accurately approximated by 
HshoAy for which we know the ground state. One finds 
that 

L 

H' = ]T [A + n (3£ - C 1 ) + {2- J V - Sl£) n k+2 ] n k 

k 



-0(£-£ 1 )^n k n k+1 (2- n k+2 ). (6) 

k 

The first term of H' can be eliminated exactly pro- 
vided that the conditions (i) A = — 0(3£ — and 
(ii) V = 2 7 0£ are satisfied. The contribution of the sec- 
ond term in eq. ^ is small since it accounts for the 
strongly diminished interaction between excited atoms 
that are at least three lattice sites apart. The third term 
vanishes exactly at £ = 1, but its contribution is negli- 
gible even away from this point since the probability for 
a simultaneous excitation of neighboring atoms is highly 
suppressed (even strictly zero in the physical subspace). 

These considerations imply that upon meeting con- 
dition (ii), i.e. for an interaction strength satisfying 
V = 2 7 0£, the ground state energy of Hamiltonian ^ 
is given by E = -£OZ where f = (1/6) [-(A/O) + 
y/l2+ (A/O) 2 ]. The latter relation is obtained directly 
from condition (i) and yields the conversion between the 
laser parameters and the square root of the fugacity, i.e. 
£. That this is indeed the case is shown in fig. [T]i where 
we compare the ground state energy Eq (red curve) with 
the numerical result (blue curve) obtained for a lattice 
with L = 20 sites. The excellent agreement indicates that 
conditions (i) and (ii) define a manifold of approximate 
Rokhsar-Kivelson points [TB] in the parameter space (O, 
A, V) were the Hamiltonian of a gas of interacting Ryd- 
berg atoms allows the approximate stochastic matrix 
form decomposition |17j shown in eq. ^ and has the 
ground state (J5J. 

We can now calculate properties of the ground state of 
the system on this manifold in the same spirit in which 
we obtained the normalization constant . Expectation 
values of classical observables such as the mean number 
of excited atoms or density-density correlations then re- 
duce to the manipulation of the partition function with 
fugacity £ 2 . The mean density of Rydberg atoms in the 
ground state is given by (N) JL = J2 k (£\n k |£) /L = 
[1 - 1/(0 + 4£ 2 )]/2 which is shown in fig. |lj). We 
can furthermore obtain the full statistics of the Rydberg 
number distribution by taking derivatives of the parti- 
tion function: The probability p k to count k Rydberg 
atoms is given by p k = [(fc!) _1 9f 2 Zj \^=o]/Z^, A common 
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FIG. 1. Comparison between the numerical results obtained 
for a lattice with L = 20 sites and 7 = 6 (blue) and the an- 
alytical expressions (dashed red), a: Energy per particle in 
the ground state, b: Mean density of Rydberg atoms on the 
lattice, c: Mandel Q-parameter of the Rydberg number distri- 
bution, d: Numerically calculated density- density correlation 
function, e: Density-density correlation function obtained for 
the state |H]). Note that at the same time as A also the po- 
tential is varied according to V = 2 7 f2£. The values of f are 
given underneath panel c. 

way for the experimental characterization of the distri- 
bution function is the Mandel Q-factor which quantifies 
the difference of the distribution p k from a Poissonian 
[2"T] . This quantity, which is plotted in fi g, [lj;, ev aluates 
to Q = {{N 2 ) - (AO 2 )/ (N) - 1 = 1/(2^/1 + 4^) _ (1 + 
8£ 2 )/(2 + 8£ 2 ). Except for £ = it is negative showing 
a pronounced sub-Poissonian behavior which is expected 
for strongly interacting systems 22 . 

A further important quantity characterizing the 
ground state is the connected density-density correla- 
tion function £i,i +m (£) = (nin 1+m ) - (n x ) (n 1+m ) = 
£7(1 + 4£ 2 )[(yT+4£2 _ 2£ 2 - l)/(2£ 2 )] m . It is shown 
in fig. [l] in panels d and e together with the numerical 
result, both again in excellent agreement. Visible correla- 
tions build up as soon as A < 0. They are exponentially 
decaying with the interparticle distance and alternating 
in sign, with anti-correlation between nearest neighbors. 
The corresponding correlation length is proportional to 
£a and reaches the system size when —A/O « 3L. 

We will now perform an analysis of the coherent prop- 
erties of the system. To this end we study the reduced 
single particle density matrix pi(£) which allows us to 
quantify the entanglement of one spin with the rest of 
the system. We find 

Pl (Z)-(l/L)^_ {N)/( L _ {N) j 

which, except for £ = 0, represent a mixed state. This 
indicates entanglement of one atom with the remaining 
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FIG. 2. a: Entanglement entropy S of a single spin with 
the remaining ones. For large negative detuning the ground 
state becomes a GHZ state and 5* reaches its maximum log 2. 
b: Density of excited atoms as a function of the detuning 
and the Rabi frequency. The black line represents the set of 
parameters where the Rydberg gas ground state is approxi- 
mately given by eq. |5}. Dashed lines are used as a guide 
to the eye delimiting the regions where the Rydberg density 
is approximately 1/3 and 1/2. c: The spin Hamiltonian Q 
can be generalized to higher dimensions (here 2d) where the 
excitation on the k-ih site (grey color) blocks the excitation 
of all sites contained in the set Gk ■ 



others which can be quantified by the entanglement en- 
tropy S = — Trpi(£) logpi(£). This function is shown in 
fig. [2^i. For large positive detuning, i.e. £ ps 0, the ground 
state is a product state |init) = J7 fc \g) k and hence no en- 
tanglement is present. S increases monotonously with 
£ and saturates at a value log 2 for £ —> oo which indi- 
cates maximal entanglement. Here the ground state is 
formally given by a GHZ state, which is the coherent su- 
perposition of the two possible anti-ferromagnetic states, 
i.e. |GHZ) = (1/V2)[\UU ■•■> + \iUt ■■■>] (even number 
of sites assumed). 

The above considerations indicate that the typical ex- 
perimental initial state |init) (no Rydberg atoms present) 
can be adiabatically connected to the fully entangled 
GHZ state by varying £ from zero to infinity, i.e. by 
varying f2 and A in time. Experimentally this is usu- 
ally done at fixed interaction strength V. The approxi- 
mate manifold of Rokhsar-Kivelson points is then given 
through {2~>n/V) 2 -(2~<A/V) = 3 which is obtained from 
(i) and (ii) and shown as the black curve in fig. ^p. The 
GHZ state is obtained by initially choosing a large pos- 
itive detuning and following this curve until one reaches 
A m i n = — 3/2 7 V, i.e. f£ m in = 0. Performing this pro- 
cess adiabatically becomes increasingly difficult as the 
number of particles increases due to an ever closing en- 
ergy gap. Eventually, this will lead to symmetry break- 
ing which singles out one of the two anti-ferromagnetic 
states or leads to domain formation. Experiments have 
to be carried out on a time shorter than the lifetime of the 
atomic Rydberg state (typically 100 /is for Rubidium and 
a principal quantum number in the range n = 40. ..70). It 



is indeed possible to find experimental parameters that 
achieve that (see Refs. pHT4l 123]). 

Let us finally discuss the generalization of Hamiltonian 
([3]) to higher dimensions and blockade ranges that can go 
beyond the nearest neighbors. To this end we replace the 
product Pk-iPk+i by an operator which projects onto 
the state n gc G fc \^)q- Here Gk is a set that contains the 
indices of lattice sites that surround the fe-th site, i.e. 
that are blocked when spin k is excited (see fig. |2b). 
The ground state of this Hamiltonian is then constructed 
analogous to the state ([5]) with the constraint being that 
a simultaneous excitation on site k and on any of the 
sites contained in Gk is forbidden. Calculations of ex- 
pectation values here again reduce to the manipulation 
of a partition sum of a classical system of hard objects. 
It is not immediately evident whether such models actu- 
ally represent an experimentally relevant system. This 
depends on whether conditions similar to (i) and (ii) can 
be found which cancel the unwanted many-body terms in 
H' . However, the knowledge of the ground state is valu- 
able, e.g. for performing perturbation theory in order to 
move away from the exactly solvable situation. 
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